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ABSTRACT 

Geometric method based on the high-order 3D Voronoi tessellation is proposed for 
identifying the single galaxies, pairs and triplets. This approach allows to select small 
galaxy groups and isolated galaxies in different environment and find the isolated 
systems. The volume-limited sample of galaxies from the SDSS DR5 spectroscopic 
survey was used. We conclude that in such small groups as pairs and triplets the 
segregation by luminosity is clearly observed: galaxies in the isolated pairs and triplets 
are on average two times more luminous than isolated galaxies. We consider the dark 
matter content in different systems. The median values of mass-to-luminosity ratio 
are 12 Mq/Lq for the isolated pairs and 44 Mq/Lq for the isolated triplets; 7 (8) 
Mq/Lq for the most compact pairs (triplets). We found also that systems in the 
denser environment have greater rms velocity and mass-to-luminosity ratio. 

Key words: Galaxies: general - galaxies: kinematics and dynamics, dark matter 



1 INTRODUCTION 

Physical properties of galaxies depend on the formation 
conditions and evolution. In addition to intrinsic evolution, 
galaxies are exposed to environmental influence (among oth- 
ers Dressier 1980, Lewis et al. 2002, Gomez et al. 2003, 
Einasto et al. 2003, Blanton et al. 2005, Weinmann et al. 
2006, Martinez et al. 2006, Park et al. 2007, 2008). The 
density of galaxies (number of galaxies per volume unit) 
or luminosity density as well as the number of galaxies in 
group/cluster or distance to the nearest galaxy is often im- 
plied as regarding the environment. The influence of an en- 
vironment can be found till 1 Mpc and even farther (Kauff- 
mann et al. 2004, Blanton et al. 2005, Park et al. 2007, 2008), 
where small galaxy groups are observed. The study of "en- 
viromental effects" in such poor galaxy groups is helpful for 
understanding the galaxy evolution on intermediate scales 
between isolated galaxies and rich groups/clusters. 

The isolated galaxies that have not sufficiently under- 
gone the influence of environment allow us to consider them 
as the " autonomous laboratories" for studying evolutionary 
processes in the galaxies. Individual properties of isolated 
galaxies (mass, luminosity, morphology, colour-index etc.) 
can be served as the standard when studying galaxies in the 
different environment (see, for example, Karachentseva 1973 
(KIG), Prada et al. 2003, Reda et al. 2004, Stocke et al. 2004, 
Verley et al. 2007) . For studing of galaxies' properties in dif- 



ferent environment it is also important to define the galaxy's 
isolation degree which is suitable to describe by some pa- 
rameter. For example, Karachentsev and Kasparova (2005) 
used the tidal index for each galaxy to study global proper- 
ties of nearest galaxies in different environments. Verley et 
al. (2007) quantified the isolation degree of KIG galaxies by 
two parameters: local number density and tidal strength. 

It is known that with the increase of galaxy systems 
richness from individual galaxies to clusters, their mass in- 
creases more quickly than luminosity (Karachentsev et al. 
1966, Girardi et al. 2002). The dark matter in small groups 
seems to be distributed in the whole volume of system in the 
case of compact groups and to be concentrated in the halo of 
individual galaxies in the case of loose groups (Mulchaey et 
al. 2003, Melnyk & Vavilova 2006, Da Rocha et al. 2008). At 
the same time the amount of dark matter in galaxy groups 
is not enough for standard cosmological model (Makarov & 
Karachentsev 2007). 

As a rule, for identifying groups by different selection 
methods, the principle of overdensity in comparison with 
background is used. The richer the group population, the 
overdensity is more strong and therefore more likely that 
such a group is physically bound. Poor groups identification 
depends strongly on the limiting parameters of the method. 
These systems can be easily confused with the random phys- 
ically not bound systems. For this reason, many authors pre- 
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fer to study rich groups/clusters only. That is why the elab- 
oration of the reliable methods for identifying small groups 
will allows us to pay more attention these structures. 

Karachentsev (1972, 1987) and Karachentseva et al. 
(1979, 2000) used 2D method for pairs and triplets selec- 
tion taking into account the isolation degree in comparison 
with the foreground and the background. The environment 
of each pair and triplet was inspected using POSSI, POSSII 
and ESO/SERC plates. This method selects effectively the 
isolated systems, which are compact in a projection with 
characteristic distance between galaxies J? ~ 50 — 100 kpc. 
With appearance of the large databases and surveys (LEDA, 
NED, SDSS, CfA2, SSRS2, DEEP2, 2dF) and information 
about the radial velocities of galaxies the application of 
three-dimensional methods of selection became generally ac- 
cepted. 3D methods of group's identification use as a rule 
two limited parameters: projected distance between galaxies 
R and radial velocity difference AV^. For example. Barton et 
al. (2000, 2003), Geller et al. (2006), Woods et al. (2006) se- 
lected close pairs of galaxies from CfA2 catalogue with R < 
50 kpc and /S.V <1000 km s~^, isolation criterion was ig- 
nored. Patton et al. (2000) investigated galaxy pairs from 
SSRS2. They considered all pairs with R <100 kpc and 
didn't find pairs with the signs of interaction with AV^ > 
600 km s^^, but most of pairs with interaction signs had 
-R < 20 kpc and Al/ < 500 km s~^ (see also Patton et al. 
2002, 2005, Xu et al. 2004, De Propris et al. 2007). Nikolic 
et al. (2004) obtained similar results for SDSS pairs. The 
authors found that a star-formation rate is significantly en- 
hanced for the projected separations less then 30 kpc. Lam- 
bas et al. (2003) and Alonso et al. (2004) considered galaxy 
pairs with i? < 1 Mpc and Al/ < 1000 km s"^ from 2dF sur- 
vey. They analyzed the star formation activity in the pairs 
as a function of both relative projected distance and relative 
radial velocity and found that the star formation activity in 
galaxy pairs is significantly enhanced over that in isolated 
galaxies at similar redshifts for R < 25 kpc and AV < 100 
km s~^. Soares (2007) showed that more than half of the 
simulated pairs with projected distance _R ^ 50 kpc have 
3D separations greater than 50 kpc. Therefore for the selec- 
tion of real physical pair it is important to take into account 
the signs of interaction, but that is not always possible. 

Triplets of galaxies are less studied than pairs. Mainly 
triplets have been selected from the catalogues of groups of 
different population (Trofimov & Chernin, 1995). For exam- 
ple, Vavilova et al. (2005), Karachentseva et al. (2005), Mel- 
nyk (2006), Melnyk & Vavilova (2006) compared the kine- 
matic, dynamical, configurational and morphological prop- 
erties of triplets from the samples formed by different meth- 
ods. They showed that physical properties of galaxy groups 
strongly depend on selection criteria. 

The main goal of our paper is to provide the uniform 
selection of the single galaxies, pairs and triplets from Sloan 
Digital Sky Survejy for later analysis of their properties: 
mass-to-luminosity ratio, colour indices, morphology and 
others. Such investigation can be helpful for study the envi- 
ronmental influence and dark matter content on small scales 
(galaxy - pair - triplet). To meet this goal we made not only 
selection of the most isolated galaxies and systems, but also 
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single galaxies, pairs and triplets with different isolation de- 
grees. 

Most authors mentioned above used simple selection al- 
ghorithms for pair identification: they consider all pairs with 
the fixed limitation parameters R and AV , which describe 
properties of galaxy pair only as a separate system, with- 
out information about their neighbours. Contrary to this we 
propose another approach. We use the second-order Voronoi 
tessellation for galaxy pairs identification where the funda- 
mental elements are pairs and the third-order Voronoi tes- 
sellation for galaxy triplets identification where the funda- 
mental elements are triplets. The geometric properties of 
the high-order tessellations give information about relative 
location of neighbouring galaxies. This allows us to anal- 
yse the correlations of group properties with environment. 
The high-order (second or third -order) Voronoi tessellation 
method has not been applied earlier for groups detection, 
unlike the first-order Voronoi tessellation. Using the galax- 
ies from SDSS DR5 survey with known redshifts allows us 
to realize the 3D approach. 

The outline of the paper is as follows. In Section 2 we 
present the method and parameters. The sample of galaxies 
is described in Section 3. In Section 4 we discuss our results 
and compare them with other works. Our conclusions are 
presented in Section 5. 



2 THE METHOD AND PARAMETERS 

Voronoi diagram was introduced by G. Voronoi in 1908. The 
most popular is the first-order Voronoi diagram (so called 
Voronoi tessellation). Fig. [1^. It is a geometric method of 
space partition on regions. Each region consists of one nu- 
cleus and all the points of space that are closer to a given 
nucleus than to other nuclei (Matsuda & Shima 1984, Lin- 
denbergh 2002). Kiang (1966) found analytic function of 
Voronoi cell volume distribution for random points in 2D 
and 3D space. A Voronoi tessellation is used widely in astro- 
physics, especially for modelling the galaxy and voids large- 
scale structure distribution (Icke & van de Weygaert 1987, 
van de Weygaert & Icke 1989, van de Weygaert 1994, Za- 
ninetti 2006, van de Weygaert & Schaap 2007), for studying 
periodicities in deep pencil-beam sky surveys (van de Wey- 
gaert 1991, Williams et al. 1991, Ikeuchi & Turner 1991, 
SubbaRao & Szalay 1992, Gonzales et al. 2000), extragalac- 
tic radiosources distribution (Benn & Wall 1995), for model- 
ing of the cosmic microwave background anisotropy (Coles & 
Barrow 1990, SubbaRao et al. 1994) etc. Ebeling & Wieden- 
mann (1993) were the first who applied Voronoi tessellation 
for finding groups and clusters of galaxies in 2D case. Later 
such approach was used by Ramella et al. (2001), Kim et al. 
(2002), Lopes et al. (2004), Barrena et al. (2005), Panko & 
Flin (2006). 3D Voronoi tessellation for galaxy groups iden- 
tification was realized by Marinoni et al. (2002) and Cooper 
et al. (2005). The application of 3D Voronoi tessellation to 
DEEP2 survey was introduced by Gerke et al. (2005) . Mel- 
nyk et al. (2006) applied the first-order 3D Voronoi tessel- 
lation for the identification of galaxy groups in the Local 
Supercluster. The authors demonstrated that the first-order 
tessellation is more useful for searching the rich clusters of 
galaxies than the small groups. In the first-order Voronoi tes- 
selation the key parameter is the volume of galaxy's Voronoi 
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cell V . This parameter characterizes a galaxy environmental 
density. The condition of cluster/group membership of cer- 
tain galaxy is the relatively small value of V . This condition 
is true when the galaxy is surrounded by close neighbouring 
galax;ies. That is why the first order Voronoi tessellation is 
not corrected for the identifing small isolated galaxy systems 
(see Melnyk et al. 2006 for details). In this paper we propose 
the high-order 3D Voronoi tessellation method for pairs and 
triplets identification. Let us consider this method below. 

Contrary to the first-order tessellation (Fig. [TJi) the 
second-order tessellation for set S distribution of nuclei is 
the partition of the space which associates a region Vi,2 with 
each pairs of nuclei 1 and 2 from S in such a way that all 
points in Vi,2 are closer to 1 and 2 than other nuclei from 
5 (Fig. [T}d). Region Vi,2 is a common cell for nuclei 1 and 
2l3 In such a way the second-order Voronoi tessellation is 
available for the identification of single galaxies and pairs. 

The third-order Voronoi tessellation is appropriate for 
the triplets identification. It is the partition of the space 
which associates a region Vi,2,3 with each triplet of nuclei 
1, 2, 3 in such a way that all points in Vi,2,3 are closer to 
nuclei 1, 2, 3 than other nuclei from S (Lindenbergh 2002), 
Fig.IB. 

Since we work with the sample of galaxies we use galax- 
ies as the nuclei of the Voronoi tessellation taking into ac- 
count equatorial coordinates a, S and radial velocities of 
galaxies Vh only. For the constructing of the 3D Voronoi 
tessellations it is necessary to determine the distances in 
3D space. Spatial distance between two galaxies one can 
split on projected (tangential) distance r and radial v (dif- 
ference of the radial velocities). We can determine the pro- 
jected distance with a relatively high accuracy, while radial 
component has errors owing to inaccuracy of radial velocity 
measurement of each galaxies and existing strong peculiar 
velocities (due to virial motions of galaxies in groups and 
clusters). As a result the galaxy distribution in space of ra- 
dial velocities is extended along radial component, so-called 
fingers-of-god efi^ect. This is attributed to the random veloc- 
ity dispersions in galaxy volume-limited sample that deviate 
a galaxy's velocity from pure Hubble fiow, stretching out a 
group of galaxies in redshift space (Jackson 1972). Various 
authors take into account this effect in their own way in 
dependence on the specificity of their problem. For exam- 
ple Marinoni et al. (2002) chose some cylindrical window of 
clusterization which extended along the radial component. 
We propose a different way that is based on introduction of 
weight for a radial component (see Appendix A). Such ap- 
proach allows us to avoid a problem of tangential and radial 
distance inequivalence and to apply high-order 3D Voronoi 
tessellation method. 



2.1 Second-order Voronoi tessellation. Pairs and 
single galaxies 

The second-order Voronoi tessellation was applied for the 
identification of pairs and single galaxies in the following 
way. Each galaxy i from set S forms the common cells with 



However it is not necessary that these nuclei lie in the common 
cell. For example nuclei 1 and 5 create the common cell Vi^s and 
they do not lie in this cell. 



the certain number of neighbouring galaxies (Fig. [TJa). So 
under neighbouring galaxies of galaxy i we understand only 
galaxies which create common cells with this galaxy. For ex- 
ample, galaxy 1 creates only 4 common cells (Vi,2, Vi,3, Vi,4, 
Vi,5) with neighbouring galaxies 2, 3, 4 and 5, respectively. 
Each pair of galaxies i,j is characterized by the dimension- 
less parameters Pij: 



Pi,j = 



(1) 



where D - space dimension, Vij - the area (for 2D) or volume 
(for 3D) of cell, m^j - distance between galaxies i and j (see 
Appendix A). 

Each galaxy has the set of Pi^j parameters. For exam- 
ple, galaxy 1 has the set pi, 2, pi, 3, pi, 4, pi, 5. We choose 
the maximum value from the set of parameters for galaxy 
1: Pmax(l) = maa;(pi,2,pi,3,Pi,4,Pi,5)- In general case for 
galaxy i: Pmaxii) = max{pij) where j assume numbers of k 
neighbours. 

We named the geometric pair in the second-order 
Voronoi tessellation such a pair which contains two galaxies 
with the common cell that have the same Pmax parameters, 
Pmax{i) ~ Pmax{2) — p. The parameter p characterizes a 
degree of geometric pair isolation. Under degree of pair iso- 
lation we understand the overdensity under background. We 
used the principle which was applied by Karachentsev (1987) 
for the selection of isolated galaxy pairs: if all galaxies have 
the same angular diameter, a neighbour of a pair of galaxies 
should be located on the sky more than five times farther 
from a pair than the separation of the pair members from 
each other. In our 3D approach the degree of pair isolation 
is described by the relation between volume Vi,2 of com- 
mon cell (it characterizes distance to neighbours) and the 
distance mi, 2 between pair members according to (1). For 
example, strongly isolated pair has a large value of p due to 
a large Vi,2 and small mi, 2 (1), see also Fig. [2^. 

We also introduce the parameter pe which describe only 
pair environment and does not depend on distance between 
pair members directly. We defined it as the mean value of 
Pj(l) and pi{2) parameters of the first and second galaxy, 
excepting p from both sets: 



Pe = 



E -=2^.(1) + EL2P'(2) 



(2) 



k + n-2 

where k and n - number of neighbouring galaxies for 1 and 2 
galaxies of geometric pair, respectively. We start sums from 
j — 2 and I — 2 for excepting 2 • p, because the first galaxy 
is neighbour for the second galaxy and vice versa. There- 
fore k + n-2 is sum of neighbouring galaxies of pair members 
excepting of pair galaxies as neighbouring for each other. 
Parameter pe depends on the distribution of neighbouring 
galaxies. Small value of pe points out that such pair is lo- 
cated in loose environment. In such case the average vol- 
ume of common cells of pair components with neighbouring 
galaxies is relatively small, and distance between them is 
large, see formula (1) and Fig. [2^. 

Single galaxy is a galaxy which is not member of any 
geometric pair. The single galaxies are field galaxies in the 
environment of geometric pairs. Every single galaxy has the 
own neighbours, single galaxies and geometric pair members 
can be among them. According to the second-order Voronoi 
tessellation the larger is the degree of galaxy isolation, the 
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greater is the number of neighbours (see Fig. [T]d in compar- 
ison with Fig. [2j3), but these neighbours locate farther. The 
best parameter which describes the isolation degree of the 
single galaxy is the mean value of all parameters pj of this 
galaxy: 



k 



(3) 



where k is the number of neighbours. Therefore the smaller 
is s value, the more isolated is the single galaxy. 



the same time, if the triplet neighbours locate far from it, 
parameter te has the small value. 

Parameters p, s, t are the basic ones and define the iso- 
lation degree of galaxy pair, single galaxy or triplet in com- 
parison with background, respectively. Parameters Pe and 
are additional ones and contain information about the 
distribution of the neighbouring galaxies (environment). 

Similarly to the second- and third-order Voronoi tessel- 
lation it is possible to apply more high-order Voronoi tessel- 
lation for the identification of galaxy quartets, quintets and 



2.2 Third-order Voronoi tessellation. Triplets 

Method of the third-order Voronoi tessellation can be intro- 
duced the same as the second-order approach (Fig. [l];). All 
points of the common triplet's cell are the most closer to 
galaxies of this triplet than to other galaxies. Similarly to 
the parameter p^j for pairs we set up the parameter 
for triplets: 



ii.j.u 



(4) 



max{mij,mi^u,mj,u) ' 

where D - space dimension, Vij^u - the area (for 2D) or 
volume (for 3D) of cell, mi,j, rrii^u, m-j.u - distances between 
galaxies in the triplet. 

Geometric triplet in the third-order Voronoi tessellation 
contains three galaxies that have the common cell and the 
same maximal parameters tmaxil) = tmax{'i) = tmaxii) = 
t. The parameter t characterizes a degree of geometric triplet 
isolation. We defined parameter of triplet environment t^ as 
the mean value of parameters ti(l), tj{2) and tu{3), except 
t from three sets: 

^ o ' (5) 

k + n + q — S 

here in case of third-order Voronoi tessellation k, n and q 
denote the number of neighbouring triplets which contain 
galaxies 1, 2 and 3, respectively. Therefore k + n + q—3 is 
number of neighbouring triplets for certain triplet, which 
contain at least one galaxy from this triplet. 

It can be seen from Fig. ^ and (4), (5) that for the 
triplet with highest degree of standing out against a back- 
ground, the isolation parameter t has the highest value. At 



3 THE SAMPLE 

For our investigation we used Northern part of the SDSS 
DR5 spectroscopic survey. Our sample is volume-limited and 
consist of objects that classified as galaxies. The primary 
sample had contained approximately 11000 galaxies with ra- 
dial velocities from 2500 km s"^ to 10000 km s~\ Ho = 75 
km s~^ Mpc~^. It is known that compliteness of SDSS is 
poor for the bright galaxies because of spectroscopic selec- 
tion criteria and the diffuculty of obtaining correct photome- 
try for object with large angular size. We tried to decrease of 
this effect's infiuence due to limiting of our sample Vh > 2500 
km s~^, i.e. in the way not to take into account nearest ob- 
jects with the large angular diameter. Such a volume limiting 
also helps us to avoid infiuence of Virgo cluster where strong 
peculiar motion exists. We checked additionally all pairs of 
galaxies with a small angular resolution and excluded identi- 
cal objects (parts of galaxies), which are presented twice and 
more in SDSS survey. All galaxy velocities Vh were corrected 
for the Local Group centroid Vlg accordingly to Karachent- 
sev & Makarov (1996). When we had applied the high-order 
Voronoi tessellation method to SDSS catalogue we limited 
our sample 3000 km s"^ < Vlg 9500 km s'^. We did not 
consider also galaxies that located within 2° near borders, 
because the correct estimation of Voronoi cell volume is not 
possible in this case. Our volume-limited sample is complete 
up to 17.7™ but contains also about 100 more fainter galax- 
ies. Final number of galaxies in the sample is 6786. 
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Figure 2. Different configurations of the galaxies: isolated close pair a) and isolated single galaxy b) in the second-order tessellation; 
isolated close triplet in the third-order tessellation c). 



Table 1. Quartiles of isolation parameters. 





Parameters 


Qi 


Q2 


Qs 


Singles 


s 


0.57 


0.77 


1.08 


Pairs 


P 


3.59 


5.72 


10.17 


Triplets 


t 


1.90 


2.62 


3.90 



4 THE RESULTS 

We applied the second-order 3D Voronoi tessellation to our 
sample 6786 galaxies and obtained 2196 geometric pairs and 
2394 single galaxies (65% galaxies of whole sample are in 
pairs and 35% are singles). We divided our samples of geo- 
metric pairs and singles on four equal parts by the parameter 
of isolation p and s, respectively. A quater of each sample 
with the highest isolation degree we called as isolated (549 
pairs and 598 singles). It means that the isolated pairs and 
singles have p > Qs and s < Qi, respectively, Qa is third 
quartile and Qi is the first quartile. See values of quartiles 
in Table [U 

Fig. |3] shows the distributions of neighbouring galaxies 
number for galaxies that are members of pairs and singles. 
We can see from Fig.[3]that number of neighbours is varried 
through range from 4 to 30 galaxies. Isolated galaxies have 
more neighbours than galaxies in other samples in average. 
It is a feature of the second-order Voronoi tessellation (see 
above) . Galaxies in isolated pairs have approximately less by 
half neighbours than isolated galaxies because of the pair's 
neighbours distribute among two members. 

Independently we applyed the third-order 3D Voronoi 
tessellation to our galaxy sample and obtained 1182 geo- 
metric triplets (52% of whole sample). The quater (297) of 
triplet sample with the highest isolation degree t > we 
called as isolated. Values of all quartiles for triplets can be 
found in Table [1] The distribution of number of neighbour- 
ing triplets is drawn in Fig. |4l As can be seen, this picture 
looks the same as for distribution of neighbouring galax- 
ies number in case of the second-order Voronoi tessellation 
(Fig.©. 




1 I I I ii 1 1 1 . I . I ' 

5 10 15 20 25 30 

Number of neigfibours 

Figure 3. Distributions of neighbouring galaxies number: 1) all 
galaxies of the sample, 2) galaxies in geometric pairs, 3) single 
galaxies, 4) galaxies in isolated pairs, 5) isolated galaxies. 



4.1 Comparison with other samples 

We cross-correlated our results with other samples. In the 
first order we compared pairs and triplets of our sample 
with Tago et al. (2008) groups selected by modified friends- 
of-friends method using the same release of SDSS. From 965 
galaxies of Tago's pairs that are located in our region of 
investigation, 686 galaxies (343 pairs or 71%) coincide with 
our pairs with difi'erent isolation degree. Median value of iso- 
lation parameter p is 9.89 that does not differ strongly from 
our isolation limit (Q3=10.17). Other 121 galaxies coincide 
with one of our pair members, median value of p for these 
pairs is 3.94. We can conclude that they locate in the field. 
The remaining 158 galaxies coincide with our single galaxies 
having isolation degree s=0.87 i.e. the also associated with 
the field galaxies (our isolation limit is less than Qi = 0.57). 
63 triplets (51%) of Tago et al. (2008) triplets that fall in our 
region coincide with our triplets. Median value of isolation 
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Figure 4. Distributions of neighbouring triplets number: 1) all 
triplets of the sample, 2) geometric triplets, 3) isolated triplets. 



Table 2. Physical properties of pairs and triplets in dependence 
on isolation parameters. 



Pairs 


N 


Sv 




Pairs 


All geometric 


2196 


99+18 


354+336 
•Ja^_i99 


V>Ql 


1677 




9c:9+217 


V>Q2 


1038 




182+if 


P>Q-i 


549 




10611? 


Triplets 


All geometric 


1182 




Qqn+333 


t > Qi 


893 


38l« 


o.,o+216 


t > Q2 


582 




97c: + 194 


t > Qa 


297 


01+11 


91 9+132 



degree for these triplets t — 4.83 is higher than our isolation 
limit g3=3.90. 

We also compared our pairs with 28 isolated Karachent- 
sev's (1987) pairs that located in our region of investigation: 
16 pairs coincide with our pairs (median value of isolation 
degree p — 21.16 is very high). Among 7 our "pairs" one 
of the component is interacting galaxy which corresponds 
to Karachentsev's pair, second component is more fainter 
galaxy, these "pairs" are not isolated (median value of p is 
3.95). Such cases appear because of difficulties in spectro- 
scopic measurement of interacting (very close) galaxies in 
SDSS. The remaining 5 pairs coincide by components with 
different pairs and single galaxies with small isolation degree 
because of these pairs are surrounded by fainter satellites. 
We found that 10 from 36 KIG galaxies coincided with our 
single galaxies, but 24 KIG galaxies fall in our pairs with dif- 
ferent isolation degree (median value is not very high p=5.8). 
Such difference can be explained by the presence of fainter 
galaxy in immediate vicinity of KIG galaxy. Actually, all 34 
isolated galaxies with fainter satellites from paper Prada et 
al. (2003), which locate in our region, coincide with our pairs 
and triplets. 

We can conclude: our singles, pairs and triplets of galax- 
ies obtained by the high-order Voronoi tessellation are in a 
good agreement with other samples, especially it concerns 
the systems with the high isolation degrees. Some lack of 
coincidence can be explained by differences in primary cata- 
logues, for example, by the magnitude depth or spectroscopy 
of interacting galaxies in SDSS. 



4.2 Isolation and main parameters 

In addition to the isolation and environmental parameters 
p, Pe, s, t and ie, we calculated also other characteristics 
(Karachentsev et al. 1989, Vavilova et al. 2005): root-mean- 
square (rms) velocity of galaxies with respect to the group 

centre (in km s'^), S„ = [^Y^L.^Vto - {VLG)f]''\ 
Sv = AV/2 for pairs, A'^ is the number of galaxies in 
the group; harmonic mean radii of the system (in kpc) - 



R 



IE 



i,k^ik ' where Rik = Xi^k{VLG}HQ^ and 
Xik is the relative angular distance; R - maximal distance 
between galaxies in triplets, for galaxy pair Rh = R; dimen- 
sionless crossing time 

T = 2HoRh/Sv (6) 
expressed in units of the cosmological time Hg^; virial mass 



3tvN{N - l)-^G~^SvRh 



(7) 



in Mq; galaxy luminosity Lr, which corresponds to Pet- 
rosian magnitude in r-band, Lq; mass-to-luminosity ratio 

Mvir/L = Mvir/J2Lr in Mq/Lq. 

First of all we investigated how the values of rms veloc- 
ity and mean harmonic radius of pairs and triplets depend 
on their isolation degree. Table [5] presents the median values 
and quartiles of Sv and Rh for pairs and triplets in depen- 
dence on isolation parameters p and t: all geometric pairs, 
pairs with p > Qi, p > Q2, p > Qa; all geometric triplets, 
triplets with t > Qi, t > Q2, t > Qa, where Qi, Q2, Qa 
is the first, second (median value), third quartiles, respec- 
tively. See values of p and t quartiles in Table [T] It can be 
seen from Table [2] that all geometric samples have largest 
values of Sv and Rh and largest interquartile range (IQR). 
With the increase of isolation degree, the medians and IQR 
of Sv and Rh decrease. Pairs and triplets selected by dynam- 
ical method (Makarov & Karachentsev 2000) have medians 



Sv 



24Zii km s" 



Rh = nOZill kpc and 41+^^ km s" 



191_8g kpc, respectively. The sample of triplets from Vav- 
ilova et al. (2005) has medians Sv ~ 30tj2 km s~^ and 
Rh = 160^77^ kpc. These values better agree with samples 
which have p,t > Q2 (i.e. p > 5.72 and t >2.62 for pairs and 
triplets, respectively). 

We also considered how the isolation degree of pairs 
and triplets depends on Sv and Rh- Tableland |4] give 
values of parameters p, pe, t, te for pairs and triplets at the 
fixed values of Sv and Rh, which correspond to intervals: 
5*1,, Rh < Qi; Qi < Sv, Rh < Q2', Q2 < Sv,Rh < Qa', 
Sv , Rh > Qa; All Sv, Rh- 

Most compact pairs with Rh less than the first quartile 
have the medians of p by 2.5 times greater and IQR almost 



Identifying Isolated galaxies, Pairs and Triplets 7 



Table 3. Parameters p and pe at fixed values Sy and Rf^ for geometric pairs. 



Pairs Sy < 10 10 < 5„ < 22 22 < S^, < 40 S„ > 40 All S„ 

P 27.6lti5:g 18.90+_l',r n.SAtltl e.utlH 14.67tl',f 

R,, < 155 pe 0.90+»:3« 0.94+0 j5 0.87+":^ 0.63+°;}:^ 0.85+^;:^^ 

TV 149 153 144 103 549 

~P S.S5+_l:Z 8-47^^:S 7.051^:^^ 7.36+^-« 

155 ^ < 354 pe 0.83l°;J^ O-SS+pJo O-S^^als 0-66^o.i2 O-'^S-a^s 

Af 146 154 133 116 549 

I 4^+1. 65 = 23+^-^* 4 70+1.79 4 (,0+1.31 4 Q4+1-S4 

354 ^ Rh < 690 pe 0.64to:i3 0.63+° ;[q 0-69lai2 0-59^o.i3 0-64^o:l3 

N 153 130 139 127 549 

o-,2+132 007+1.21 onc + 1 17 2 42+°-^° 2S3+1-10 

p o.lz_g g2 •J-^'_0.81 '>-"'J-0.70 ^-^^-0.43 ^•"'-'-0.65 

i?.,^690 p. 0A7tl:ll OAOt^ol 0.521^^ 0.47«:J« 0.481°-° 

N 101 112 133 203 549 

~P 7.27+_rrt 7.57t«3:S 5.63l?;^I 3.8ltl:g 5.72t^:« 

Alii?, p, 0.7ll°:25 o.70t°;g 0.69t°:f, 0.55«;}^ 0.651°;?^ 

A 549 549 549 549 2196 



Table 4. Parameters t and te at fixed values Si, and for geometric triplets. 



Triplets 5„ < 26 26 < S„ < 41 41< 5„ < 61 5„ > 61 All Sv 

, 92+2.24 . 00+2. 84 n 47+I.4O , 9c:+0.71 o 7^+2.53 

I O.ZZ_2 j2 ^•''■^-1.42 ■^■^'-1.11 ^■^''-0.29 '^•'■^-1.23 

Rh < 205 te 0.43_o 0.41_Q pg 0.36_q gg 0.29_g gg 0.38_q g^ 

N 91 92 70 43 296 

, 044+1.74 o OR+I.IS 9 C4+1-05 9 2c:+0.69 01 -,+1.17 

t o.iy:_g 32 ■-'•■-'O_i.04 ^■"^-1.64 ^■^•-'-0.29 '-'•^^-0.85 

205 ^i?.„< 390 ie 0.41+°:05 0.38+°:°5 0.37+°'°* 0.32+0;°3 0.37+°:06 

N 81 87 83 44 295 

. 2 Qfi+"'83 9 47+0.92 , ^n+O-^" 9 9 49+0-79 

390 < 723 0.3lt°:0° 0.32l0:g5 0.33t° °6 0.29t° °^ 0.32t°:°« 

N 71 65 76 83 295 

, 1 09+0.92 , 07+0.33 1 SQ+0.61 1 74+O.69 1 SH+0-64 

i 3g J^-^'_0.25 ^■°^-0.35 ^-'^-0.31 ^•"°-0.38 

> 723 to 2^+°-°'^ 27+°-°3 97+0-04 n 94+0-02 n 2'^+°-°^ 

rSh ^ ii-i te U.ZO_g g5 U.Z(_g g5 U.^r_g Q4 U.^4_g 03 U.ZO_g g4 

N 53 51 66 126 296 

, 090 + 1.96 •5nS+l-37 9 f;t;+1.02 9 no+0.57 9^.9+1.28 

C (j.^ J.Ue5_j^gj^ Z.DO_g gg Z . _ g 44 Z.DZ_g>,-^ 

All D * r, 07+0.09 n •^t;+0'08 n •^'?+0-07 ^ 97+0-05 r. 09+0.07 

All Ry, ie 0.37_g gg 0.35_q g^ 0.33_g gg 0.27_g g4 0.3i_g gg 

N 296 295 295 296 1182 



by three times greater than the values for sample "All Sv, 
Rh" ■ That means that compact pairs are more isolated (or 
more contrast ones in comparison with the background) than 
in average galaxies in the whole sample, but also they are 
characterized by different isolation degree. At Rh < Qi and 
p,t> Q3, 65 % pairs (47 % triplets) coincide. 

Parameter pe for the compact pairs has a large disper- 
sion also. This means that similar compact pairs are both 
in densest and loose environment. Thus, the isolation degree 
and IQR change depending on the values of rms velocities: 
the most compact pairs with rms velocity less than the first 
quartile are most isolated. For these pairs parameter p al- 
most by 5 times greater than median of "All Rh, Sv" sam- 
ple and IQR by 4.5 times greater. Most compact pairs with 
Sv > Q3 have isolation degree the same as galaxies of whole 
sample. The values of isolation parameter decrease with in- 



creasing of Rh and Sv, thus, at Rh greater of third quartile 
the p becomes almost by 2 times less than the median of a 
whole sample. It means that these pairs are less isolated. It is 
necessary to note, that p significantly changes in dependence 
on Rh, than Sv But all of pairs with Sv > 40 km s~^ (Qs) 
are characterized by small isolation degree, than pairs with 
Sv < 40 km s~^. Additional parameter pe also decreases 
with Rh increasing and for Rh > Q3 it has the minimum 
value with a small variation {p is also minimal). So, all wide 
pairs {Rh > Q2) are in low dense environment. It is obvious 
that for the widest pairs to be isolated there are not enough 
free space (without galaxies) . Most probably these pairs are 
the accidental ones in the common field. From Table 3] it is 
easy to be convinced that all above mentioned tendencies 
are similar for triplets but with their own values. 
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Figure 5. Luminosity distributions of isolated galaxies (dotted 
line), galaxies in isolated pairs (solid line), galaxies in isolated 
triplets (thick line). 



Table 5. Comparison between the luminosities of the galaxies 
from different samples. 



Sample of galaxies X IO'-'(Lq) 

Median Mean (SD) 



All single 


2.23 


+3.66 
-1.15 


5.17 


(8.00) 


All geometric pairs 


2.38 


+4.27 
-1.29 


5.79 


(8.66) 


All geometric triplets 


2.48 


+1.36 
-9.14 


6.00 


(8.77) 


Isolated 


1.80 


+3.10 
-0.77 


3.52 


(4.47) 


Isolated pairs 


2.89 


+7.20 
-1.73 


6.74 


(9.46) 


Isolated triplets 


2.68 


+6.80 
-1.52 


6.53 


(9.97) 



4.3 The luminosity content 

As considering the luminosity-density relation (for example 
Park et al. 2007, 2008 and references therein) we can note 
the following: galaxies are more luminous in the high density 
regions than in the field. Karachentseva et al. (2005) showed 
that galaxies even in such poor groups as the triplets are 
more luminous than the isolated galaxies. We compared the 
luminosities of single/isolated galaxies with the luminosities 
of galaxies in pairs and triplets. Their medians with quartiles 
and mean values with standard deviations are presented in 
TableO Fig.[5]shows the distribution of galaxy luminosities. 

You can see from Table [5] and Fig. [5] that the median 
and mean values of luminosities are greater for galaxies in 
pairs and triplets than for single/isolated galaxies. The mean 
values of luminosities of isolated galaxies and galaxies in 
isolated pairs differ significantly (with the probability of > 
0.99 by t- criterium), but differences for "All" samples are 
not statistically significant. The same conclusions we made 
from comparison of the sample of single/isolated galaxies 
with galaxies in triplets. 

So, galaxies in isolated pairs and triplets are two times 
more luminous than isolated galaxies. It is necessary to note 
that our method of group identification doesn't take into ac- 
count individual physical characteristics of galaxies, there- 
fore the effects of selection are absent here. The fact that 
for " AH" samples the mean values have small difference may 
serve as the evidence of influence of non-physical (acciden- 
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Figure 6. Dependence between slope a and R: for the pairs 
log{p) = a ■ log{Mi,ir / L) + b (thin line) and for the triplets 
log(t) = a ■ log(Myir / L) + b (thick line). 
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Figure 7. Dependence between correlation coefficient in relations 
from Fig.[6]on R: for the pairs (thin line) and for the triplets (thick 
line). 

tal) groups in these samples, the wide systems are dominated 
among them. 

4.4 Mass-to-luminosity ratio 

For studying dark matter presence and distribution in 
small galaxy groups we used the mass-to-luminosity ratio 
(Mvir/L) as the quantitative indicator of dark matter con- 
tribution. We checked the mass-to-luminosity ratio in depen- 
dence on isolation degree of pairs and triplets. We plotted 
the dependences of system isolation p, t on M^ir/L in the 
narrow bins of R, because the isolation parameters p, t and 
virial mass Mvir depend on the projected distance between 
the galaxies R (see (1), (4) and (7)). Fig. [6] presents the de- 
pendence of slope a in p, t-M^ir/L relation on R, and Fig. [7] 
presents dependence of correlation coefficient on R (details 
in signatures to the figures). 

We see from Fig. [5] and 7 that if the system is located 
in denser environment, it has greater value of M^ir/ L (a < 
-0.1, correlation coefficient < -0.3). Moreover, such depen- 
dences are true only on the scales approximately up to 150- 
200 kpc for pairs and up to 250-300 kpc for triplets. For 
greater values of R the p,t — Mvir/ L dependences were not 
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observed. This conclusion only confirms our results (subsec- 
tions 4.2 and 4.3) where we showed that the widest systems 
most probably are accidental formations, because they don't 
stand out against a background. 

Our results agree with the work by Park et al. (2008) 
where authors studied the dependence of morphology and 
luminosity on environment. The strong dependence of mor- 
phology on the nearest neighbour was found at the distances 
of about 200 kpc. On the other hand authors (see for 
example Lambas et al. 2003, Alonso et al. 2004, Woods et 
al. 2006, Patton et al. 2000, De Propris et al. 2007) on the 
base of analysis of physical properties such as starformation 
rate and colour indices concluded that differences between 
galaxies in pairs and fields galaxies are significant in pairs 
with R <20-50 kpc. Dependence on velocity dispersion is 
weaker, but Lambas et al. (2003) and Alonso et al. (2004) 
found significant differences in AV < 100 km s~^. 

Because of maximal distance J? is a constant, enhance- 
ment of Mvir/L value is a result of increase of rms velocities 
in groups. In other words, if the relatively compact pairs 
and triplets located in denser environment, then they have 
greater value of mass-to-luminosity ratio due to larger virial 
motions. The same result was obtained by Einasto et al. 
(2003). The authors found that loose groups in the neigh- 
bourhood of a rich cluster are typically 2.5 times more mas- 
sive and 1.6 times more luminous than groups on average, 
and these loose groups have velocity dispersions 1.3 times 
greater than groups on average. 

Table [6] presents the medians and quartiles of system 
crossing time r and mass-to- luminosity ratio M^ir/L for 
samples of pairs and triplets that characterized by different 
isolation degree and compactness. 

From Table |6] follows that values of M^ir/L at i? < 50 
kpc for pairs and R < 100 kpc for triplets is 7-8 Mq/Lq 
and increase with enhancing of Rh faster for triplets than 
for pairs. 

The values of M^ir/L are in a good agreement with 
the corresponding values within IQR for groups which were 
selected by dynamical method. The values of M^ir/L are 
18t?4 Mq/Lq and 32t^2 Mq/Lq for pairs and triplets, re- 
spectively (Makarov & Karachentsev 2000) and for triplets 
from Vavilova et al. (2005) 35j:^!^ Mq/Lq. The median 
of most compact triplets is in agreement with median 
Mvir/L = 131]^^ Mq/Lq obtained for interacting triplets, 
where star formation rate is high (Melnyk 2006). The me- 
dian of isolated triplets is in a good agreement with median 
Mvir/L = 47^33 Mq/Lq which corresponds to the sample 
of isolated Northern (Karachentsev et al. 1979) and South- 
ern (Karachentseva & Karachentsev 2000) triplets. 

An observable agreement with previous results is addi- 
tional evidence of the correctness of applying the high-order 
Voronoi tessellation method for identification of pairs and 
triplets. 



5 CONCLUSIONS 

For the first time we introduced and applied the high-order 
3D Voronoi tessellation method for identification of isolated 
galaxies, pairs and triplets. 

We used volume-limited sample of 6786 galaxies from 
Northern part of the SDSS DR5 spectroscopic survey (3000 



km s"^ ^ Vlg ^ 9500 km s"^). We select single galaxies 
and pairs by the second-order Voronoi tessellation, as well 
as triplets by the third-order Voronoi tessellation method. 
As a result we formed 2196 geometric pairs, 1182 triplets 
and 2394 single galaxies. Then we introduced parameters p, 
t and s to characterize the isolation degree of pairs, triplets 
and single galaxies, respectively. We did not do a clear di- 
vision between physical gravitationally bound systems and 
non-physical ones, following supposition, that than system 
is more isolated, the more probability that it is physical. We 
named the quaters of our single galaxies, pairs and triplets 
with s < Qi, p > Qs, t > Q2. , respectively, as "isolated" 
(i.e. with highest isolation degree). The values of quartiles 
see in Table [1] So, we consider the properties of single galax- 
ies, pairs and triplets in dependence of their isolation degree 
(in different environment). 
Our main conclusions. 

Compact pairs {Rh < 150 kpc) and triplets {Rh < 200 
kpc) are more isolated in average than systems in geometric 
samples, thus they are characterized by different isolation 
degree. The wider pair (triplet), the smaller isolation degree 
is observed. Small values of parameters pe and te are the ev- 
idence of loose environment of these systems (they have not 
a "free space" to be the isolated groups). Thus, we consider 
wide geometric pairs and triplets as accidental ones in the 
common field. 

We compared the luminosities of single galaxies and 
galaxies in geometric pairs and triplets. It was shown that 
galaxies in isolated pairs and triplets are two times more lu- 
minous than isolated galaxies. On one hand it is the evidence 
of our geometric method accuracy. On the other hand, we 
can conclude that in such small groups as pairs and triplets 
the luminosity-density relation is observed. 

We considered dark matter content in our groups. The 
median values of M^ir/L for our samples limited by differ- 
ent criteria are 12 Mq/Lq for isolated pairs, 44 Mq/Lq 
for isolated triplets, 7 (8) Mq / Lq for most compact pairs 
(triplets) with _R < 50 (100) kpc. Note, that for most com- 
pact (close or interacting) systems is not very large differ- 
ence in dark matter content for pairs and triplets, but for 
isolated triplets the M^ir/L is larger in 3 times than for 
pairs. These results are in agreement with works of other 
authors. We also found that the pair/triplet is less isolated 
system (in denser environment), when M^ir/L greater. This 
relation testifies that galaxy systems in denser environment 
have greater rms velocity (because of M^ir ~ at fixed 
distance between galaxies). The p,t- M^ir/L dependences 
are observed only for compact systems (up to 150-200 kpc 
for pairs and up to 250-300 kpc for triplets). 

We concluded that 3D Voronoi high-order tessellation 
method is the effective tool for small groups identification 
and studying their properties in dependence on environ- 
ment. In the next paper we will present results on morpho- 
logical content and colour indices of galaxies in pairs and 
triplets in comparison with isolated galaxies and their envi- 
ronment . 
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Table 6. Medians and quartiles of the system crossing time and mass-to-luminosity ratio for pairs and triplets. 



Pairs 


N 


T 




Triplets* 


TV 


r 


Myir/L 


P>Q3 
R < 50 
R < 100 
R < 150 
R < 200 


519 
133 
335 
549 
684 


1 12+1'2* 
^■^^-0.59 

0-18_o 07 
40+" *'' 

U-^U_g 21 

u-or_Q 32 

70+°-''^ 

'J-"J_0.34 


-■2+37 

7+28 
'-6 

Q + 31 

^-8 
14+44 

17+50 
i'_14 


t > Q3 
R < 100 
R < 200 
R < 300 
R < 400 


297 
16 
70 
136 
240 


1 06+"**** 
J^-UD_(j 42 

43+f 
'^•^•^-0.24 
n 7n+0.33 
'^•"J-0.28 
n 71+0.39 
"•'^-0.30 


8l« 

20+**^ 
^'-'-13 
24+41 

01+49 
■^^-18 



* Maximal distance R correlates with Rf^: R = 1.41 ■ Rf^ + 266, the correlation coefficient is 0.87 for triplets, and R = Rf^ for pairs. 
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APPENDIX: 
DISTANCE 



RADIAL AND PROJECTED 



To avoid of fingers-of-god effect we used the correction fac- 
tor in calculation of space distance between galaxies. If the 
projected r and radial v distances were an equivalent, we 
should calculate a space distance between two galaxies as 



,2 2,2 

a = V + r . 



(8) 



Obviously we work in the space of radial velocities and un- 
der r we mean r-Ho- We inserted the certain factor A; < 1 as 
a weight of the radial component. This factor is responsible 
for the relative virial motion of galaxies. In that case the 
modified distance is = ii^fc^ + r^. After some transfor- 
mations we obtained the equation of ellipse: 



Table 7. Physical properties of pairs in dependence on isolation 
parameter and Vp. 



Pairs 


N 


Sv 


Rh 


Vp 


=100 km s^i 




All geometric 


2211 


^g+15 


oqQ+360 


p > 2.93 


1658 


17+11 

l'~9 


308^226 


p > 4.72 


1005 


14+8 


999+163 
^^^-113 


p > 8.15 


553 


12+'' 
r^„5 


149+92 


Vp 


=500 km s-i 




All geometric 


2210 


24+21 
^^-13 


QQQ + 331 

33»_i9i 


p > 3.83 


1657 




907+200 
^■^'-120 


p > 6.15 


1105 




167+"^ 

-84 


p > 11.19 


553 


-17+10 
l'-9 





(10) also changes to (8). So in these cases the virial veloci- 
ties have a week action on the distance measurement accu- 



racy. 



If d 



than 1 + 



> 1, i.e. the weight of radial 

2 



distance decreases. In case d << Vp denominator (l + -j- 
tends to infinity and the radial component loses significance, 
while projected distance r keeps it. For the our calculations 
we used Vp = 300 km s"'^ as the value of typical relative 
velocities in small galaxy systems (Ceccarelli et al. 2005). 
Moreover Karachentsev et al. (1989) showed that major- 
ity of physically bound triplets have rms velocity Sv < 300 
km s^^. 

For the testing our method robustness we compared 
main parameters in cases of different Vp. Table [7] presents 
the medians and quartiles of 5*1, and Rh for the obtained 
galaxy pairs at Vp = 100 km s~^ and Vp — 500 km s~^, re- 
spectively. The dependences Sv , Rh on Vp are not significant 
(see Table 2 for comparison). Increasing Vp by 5 times pro- 
duces the changing of rms velocity and projected distance 
between galaxies by 1.4 times only. In our paper we choose 
the compromise value Vp — 300 km s"'^. So, such approach 
let us to avoid the mistaken identification of wide pairs with 
a small velocity difference. 



-f r 



(9) 



where m is the minor semiaxis and ^ is the major semiaxis 
of the ellipse. We took into account some tolerance Vp in 
measuring v. For that we labelled a major semiaxis as ^ = 
d + Vp, where is a distance in the space of radial velocities 
according to (8). So the weight of the radial component is 
k — , . Since new distance m is a part of formula for 
k, therefore for the simplification we used here m — d and 
obtained new distance 



(10) 



In case of Vp — the formula (10) changes to (8). If the 
galaxies are located at the great distance d » Vp, formula 



